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Abstract 

During the last decade Levy processes with jumps have received increasing popularity for 
modelling market behaviour for both derviative pricing and risk management purposes. ? in- 
troduced the use of empirical likelihood methods to estimate the parameters of various diffusion 
processes via their characteristic functions which are readily avaiable in most cases. Return 
series from the market are used for estimation. In addition to the return series, there are many 
derivatives actively traded in the market whose prices also contain information about param- 
eters of the underlying process. This observation motivates us, in this paper, to combine the 
return series and the associated derivative prices observed at the market so as to provide a more 
refletive estimation with respect to the market movement and achieve a gain of efficiency. The 
usual asymptotic properties, including consistency and asymptotic normality, are established 
under suitable regularity conditions. Simulation and case studies are performed to demonstrate 
the feasibility and effectiveness of the proposed method. 

KEYWORDS: Empirical likelihood; Levy processes; Diffusion processes; Characteristic Func- 
tions. 

1 Introduction 

Brownian motion and normal distribution have been widely used in the BlackScholes option-pricing 
framework to model the return of assets. Stylised facts, however, contradict with the model assump- 
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tions specified in the Black-Scholes framework. This motivated studies to modify the Black-Scholes 
model to explain the empirical phenomena. One direction of extension is to model the asset return 
dynamics through Levy processes which are able to capture jumps and the asymmetric leptokurtic 
features. Readers are referred to ? for a full account of Levy processes and their applications in 
finance. 

Empirical likelihood, introduced by ?, provides an alternative, nonparametric approach to in- 
ference. By placing a probability pj on the jth observation and computing a profile likelihood, 
the method can be used to construct nonparametric point estimation as well as confidence regions 
for the parameters of interest. ? first linked estimating equations with empirical likelihood. In 
the paper, they developed methods of combining information about parameters in over-constrained 
optimisation problems which are frequently discussed in econometrics or finacial literature. When 
the number of estimating equations is larger than the number of parameters of interest, the empir- 
ical likelihood estimation procedure will automatically combine the constaints by assigning them 
appropriate weights and produce an efficient estimate. ? gives a good review on the development 
and applications of empirical likelihood. 

Statistical inference based on the characteristic functions was proposed by ?, ? for independent 
observations. ? suggested using characteristic functions as constraints for the empirical likelihood 
estimation. Such an approach makes use of the advantage that the characteristic functions of many 
diffusion processes are readily available, but it does not incorporate information from the market 
including derivative prices, for instance, which provide informative and most-updated knowledge 
of the parameters of interest. The key goal of this paper is to discuss how one can make use of 
the market data in the empirical likelihood estimation procedure to obtain more accurate estimates. 

Let {St}t>o be a continuous-time Levy process that records the evolution of a financial security 
over a period of time. Assuming that <S"s are observed over a collection of discrete time points: 
0,S,2d, . . . ,n5, over a time span [0,n(5], we can treat the difference of any two consecutive obser- 
vations, i.e. the increments, as a set of independent observations with the same distribution since 
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increments of a Levy process are independently and identically distributed. In other words 

Rj := log Sjs - log5'(j_i)5 ^iid Fe, say 

whose characteristic function is given by 

(p{t;e) = E^[exp{itRj};e] = J exp{itr}F0{dr), 

where P denotes the expectation taken under the physical measure and denotes the parameters 
of interest that governs the process {log St}t>o- Of course, using the maximum likelihood approach 
can produce the most efficient parameter estimates. This is, however, only possible when the den- 
sity function is readily available, which is not the case for most of the Levy processes. In this 
paper, we follow ? to formulate an estimation procedure using the empirical likelihood with char- 
acteristic functions as one of the constraints. Observe that a characteristic function contains the 
same amount of model information as what a probability density function can carry, it is sensible 
to incorporate them as one of the estimating equations. Instead of having the return sequence as 
the only source of data, wc can, in fact, incorporate information from actively-traded derivatives in 
order to provide a more timely estimate of the model parameters. In this paper, prices of European 
call options on the same underlying asset are used as moment constraints for empirical likelihood 
estimation procedure. Due to put-call parity between Euorpean calls and their put counterparts, 
it suffices to include just call prices as the put counterparts should contain the same amount of 
information. 

The remainder of the paper is organized as follows: we first define the notation and describe 
the methodology needed in Section 2. Sections 3 provides readers with specific examples on how to 
apply the results in Section 2 to carry out the estimation procedure. Section 4 extends the model 
to multiperiod case. A simulation study and a case study are given in Section 5 and 6 respectively, 
followed by a discussion in Section 7. Proofs are relegated to the Appendix. 
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2 Methodology 



2.1 Single Period Model 

Throughout this section, we assume that Ri, . . . , Rn are iid random variables with distribution F 
and the characteristic function (f){t;9), whose closed-form fomulation can be readily obtained. To 
begin, we start with the simpliest possible: In addition to the return series, we also observe risk- 
free rate r as well as a call option with maturity 5 and strike K. Levy processes have independent 
stationary increments and so the above set up fits Levy process. 

Following ? and ?, we study the maximum empirical likelihood estimator (MELE) based on 
constraints due to both the characteristic function as well as option prices as follows. First, it is 
easy to see that the equation <j){t\6) = E[e^*^i\ provides us with two constraints on G: 

n n 

y^Pj cos{tRj) = Re[(j){t; 6)] and ^^Pj sin{tRj) = Im[(j){t; 6)], 

where Re{z) and Im{z) denote respectively the real and the imaginary parts of z. 

For the option constraint, denote c{SnS,K,r,S,0) the call price observed at time nS, with the 
underlying asset price SnS and strike K that matures at (n+ 1)6. To simplify the notation, for the 
rest of the paper, we suppress the subscript 6 and use Sn and Rn to denote the underlying asset 
price and the associated return at time nS respectively. 



Observe that 



= [g-r5 max{5„+i - K, 0} - K, r, 6, 9) 
= E^-"-^ max{S„e-^+i - K, 0}|S„] - c{Sn, K, r, S, B) 
= E^ e-'^ max{5„e^"+i - K, 0}^ 0) - Z{Sn, K, r, 6, 9) 



Sn 



(1) 



where -^{R;d) represents the Radon-Nikodym derivative (or the density ratio) of R that adjusts 
the difference between the probabilities defined under the physical and a risk neutral measures. 
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Option prices are usually specified through moneyness which is denoted by m = S^jK. ([T]) can 
be rewitten as 



= E"" 
= 



e-''^ max{e^ - m, 0}^(-R; 0) - c(m, r, 6, 6) 



c{Sn,K, r, 6,9) 



Sn 



(2) 



which gives an additional constraint. Note that c(m, r, 6, 0) = c{Sn,K, r, 5, 0)/Sn is independent of 



Sn in most cases; see Section 3. 



The above derivation differs from ?'s canonical approach in which case historical returns are 
used to construct n possible values for the asset price one period from now, i.e. 

Sn+l = SnB^^J = 1, . . . ,n. 

That is, the previous realised returns are used to construct possible prices at (n + 1)5. ? used 
a similar approach to investigate alternative tilts for non-parametric pricing. They proposed the 
following estimating equation: 



E"^ 



e-''^max{Sne^ - K,0} = c{Sn,K,r,5,e 



which will, however, create bias in cases with small sample sizes because of the projected asset 
price. The difference between the magnitudes of Sn, thus the additional constraint, with the two 
constraints derived from considering the characteristic function of the return series may produce 
unstable numerical estimates. 
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2.2 Model Setup 

Denote pi{t), . . . ,Pnit) be probability weights allocated to the residuals ^)}j=i,...,n) where 



cos{tRj) -(l)^{t;e) 
sm{tRj) - (l)\t]e) 
y e-"^^ max{e^J - K,0}^{R j -,6) - c{m,r, 6, 6) J 



(3) 



An empirical likelihood for at t is given by 



Ln{r,e) = l[pj{t), 



(4) 



subject to constraints Yl^=iPjif) — 1 Yl^=iVj{t)gj{t;0) = 0. Applying Lagrange-multiplier 
approach as we usually see in maximum empirical likelihood derivation, we see that (jH) is maximized 
when 



1 



1 



n\{t-eygj{t-ey 

where A(t; 6) is a Lagrange multiplier in M'^ satisfying 



Q 



In 



9j{t;0) 



n 



Hence, the local log empirical likelihood ratio becomes 

n 

£n{t; 0) = 2 ^ log {1 + X{t- Oygjit; 6)} . 

Like ?, we consider integrating init', 0) against a probability weight 7r(t), an integrated empirical 
likelihood ratio for 6 is given by 



(6>)= / £n{t;e)TT{t)dt. 

JteM. 



The maximum empirical likelihood estimator (MELE) for 9 is defined as 



dn = argmin£„(0). 

06© 
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Remark: argmin is considered because —2 has been multiplied to the EL ratio ini^)- 



The above estimation procedure can be easily extended to situation in which there is more than 
one option traded in the market. In other words, options with the same maturity but different 
strikes can be added as additional constraints. For multiple strike constraints, say there are k 
European calls with moneynesses nii (i = 1, . . . ,k) respectively, one can simply rewrite ([3]) as 



/ cos{tRj) - <j)^{t;e) 

sm{tRj) - (p^it-e) 



\ 



\ max{e^^ - mm,0}§{Rj;e) - c{mm,r, 6, 6) ) 
through which we can obtain Qel using the same estimation procedure. 



2.3 Multiple-period Model 

The above framework can be further extended to incorporate options with different strikes as well 
as different maturities. Similar to the single-period case set-up, suppose we have observed a series 
of returns {Rj\j=\ y^ with the current asset price Sn- In addition, we can also obtain prices for 
calls with different maturities M and moneynesses m. 



The procedure will follow closely to the methodology proposed in Section 2.2. We start from 
the simplest case in which there are two groups of calls: one group contains A^i calls with different 
moneynesses but the same maturity b while the other group containing N2 calls with different 
moneynesses but the same maturity lb. 



The single-period case can be dealt as what we have done in Section 2. For the double-period 
model, we can view each pair of consecutive (non-overlapping) returns as a single observation. In 
this case, the double-period model can be reduced to the single-period model with n/2 number of 
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observations. Essentially, it means 



Rf^ =R2j-i + R2j, i = l,...,n/2 



The corresponding set of estimating equations for the calls that mature in 26 can be written as 



cos{tRf)-(l)^{t;e) 

?(2)^ 



sm{tRy')-(t>^{t;e) 
e-"-^ ms^{e^f - mf ^ , 0} § ; 0) - c(2) (mf ^ , r, §, 6) 



where c^^\ and m^^) denote respectively the call prices and their corresponding moneynesses. 



(2) (2) (2) 

We can define (t), Ql^ and In {t\0) accordingly for double-period model. In general, we 
can also extend above extension to multiple-period case in which g becomes 



max{e 



cos{tRf^)-(l)^{t;e) 
sm{tRf)-(l)^{t;e) 
■<\0}f (i2f ;0)-cW(mf),r,<5,0) 



V e-^max{e^f -i^(%,0}f (i?f ;0)-c('=)(m(%,r,<^,0) J 



where Rf^ = + ... + Rkj. 



Following the idea of ?, we try to express our overall likelihood as a sum of all the sub-empirical 
likelihood. The maximum likelihood estimator for 6 can be defined similarly as 



no 



0„ = argminV4^)(0), 



flee 



(5) 



where no denotes the number of unique maturities of the options observed. Readers should be 
noted that, for simplicity, we just use the call prices as constraints. One can use any other option 



prices as long as they can write down the estimating equations. The inclusion of the puts may not 
help estimation due to the put-call parity. This methodology, of course, performs worse when the 
maximum maturity becomes long that leads to a huge reduction of the number of observations. One 
should note that, however, only options with short maturities are traded actively. These options, 
meanwhile, provide the most up-to-date, thus useful, information about the parameters. 

3 Models and Examples 

In this section, three commonly used models with known characteristic functions are considered. 
Discretely observed data are used to investigate the performance of the proposed empirical likeli- 
hood estimator to provide an accurate estimate of the unknown parameters of the continuous time 
models studied. 

3.1 Black-Scholes Model 

Suppose the stock price Sf follow the geometric Brownian motion 



where Wt is a P-Brownian motion. Again, we denote the historical returns of the previous n 




(6) 



trading period as Rj = log Sj — logS'(j_i), j 



= 1, . . . , n, we know that for each j, 




The characteristic function of Sj is given by 



(t){t; e) = exp{Siit{n - - aH'^/2)}, 



where 6 = (/x, a). Hence, for any j = 1, . . . ,n, 



(7) 



is an estimating equation for 0. 
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In addition to the return series, we also observe option prices traded at time n5, each of them 
expires in the next period of length 6: {c(mj , r, 0)}j=i^...^fc. From these k option prices, we can 
write down an estimating equation for the parameters = {fJ-,cr): 



e-''^ max{e^ - m, 0}-^(-R; 0) - c{m, r, S, 6) 



where, if (El) holds. 



{R-e) 



1 



exp 



{R-{r- a'^/2)5f 



exp 



exp 



r — /i 



R 



This leads to the following estimating equation 



E 



max{e^ — m, 0} — c(m, r, 5, 6) | exp < — —R — ^ ~^ 



r" — ^ ^ (r — 



0. (8) 



3.2 Black-Scholes Model with Merton Jumps (BS-MJ) 



Empirical studies suggest that log return sequences usually exhibit skewness and an excess kurtosis 
(compared with a normal distribution). In order to devise a model that can provide a better fit 
to the financial market data, ?, believing that the Black-Scholes solution is not valid as the stock 
prices dynamics should not be presented by a stochastic process with a continuous path, proposed 
Black-Scholes Model with jumps (BS-MJ), which is specified as follows: 



dSt = (// - XK)Stdt + aStdWt + (Jt - l)StdNt, (9) 

where Nt is a Poisson process with intensity parameter A > and Jt is the jump size fol- 
lowing a lognormal distribution log— J\f{fi j, a j) and is independent of Wt- Xk := XE[Jt — 1] = 
A(exp{/i J + (Tj/2} — 1) is the compensator of the compound Poisson process ( Jt — l)StdNt. 
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By Ito's lemma for jump diffusion processes (see ?), ([9]) can be rewritten as 



dlog St = {fj.- Xk- (j'^/2)dt + adWt + log JtdNt, (10) 

under the physical measure P. Despite the fact that there is no closed form density for log St, its 
characteristic function is given as follows: 

</.(t; 61) = exp {5 ^it{fi -Xk- a^l) - + X (^g^'^./*— j*V2 - l)] } . (11) 

By constructing a hedging portfolio, ? proposed that the European call option price on an equity 
that follows the dynamics given by (fTOj) V{St,t) should be the solution of 

ay _.2 02 n2T/ qt^ ay 

'm^^ds^^ ''^'W mv{JtSt,t) - v{St, t)] - xSt^EiJt - 1] = 0, 

ft t 

which is equal to 

V{t,St) = e-''^E^^'[max{S(^t + 6) - K,0}\St] 

with 

A = A(1 + k) = Aexp{/xj + cr5/2} 



(Tn = \/a'^ + naj/6 



J 

I 1,1 1. I — T— I l.t 

r — Ak + 



n/ij + naj/2 



Again, we need to compute the Radon-Nikodym derivative between the two measures P and 
Qj\/. Using the inverse Fourier transform formula, we can express the density of the Merton's jump 
diffusion model under physical measure P as follows: 
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1 f °° 

fix) = — / 0(t;^)e-*-'di 



oo 
oo 



„A<5 

2vr J-oc 
27r 



gA5 



27r 

n>0 



n>0 

^ -^^ — / exp{i(5t(;u - Ak - crV2) - - itx + intfi j - nt^a'^j /2}dt. 

„-,n ^-OO 



Using the identity that e-'^^'+^^dz = e"^"/^", putting a = 5a'^ /2 + no] 12, b = 
cr^/2 — Xk)6 — X + nfij, we can write 



Kfi-cr^ /2-\K)S + ntJ.j-x]^ 



/(x) = N 



which is a fast converging sequence. So, the Radon-Nikodym derivative required is 



dF 



E 

n>0 



(A(5)" e 



2((5(t2+,ict^) 



6a^ + 



[(t — (T^/2-AK)<5+nfij-ii 



n>0 



n! 



In other words, the corresponding estimating equation that is derived from an option is given 



by 



E 



e""^ max{e^ - m, 0} - ^ 



n>0 



n! 



E 

n.>0 



(A5)" e 



[(M-tT^/2-Afi)5 + nMj--Rj]^ \ 
— — T- 



2(<5o-2 + nCT^) 



n! 



5cj2 + 



E 

n>0 



(A(^)" e 



[(r-CT2/2-Are)5 + nMj--Rj]^ \ 



-1- 



ni 



ncr' 



To generate log St from pop , we use a sequence of Bernouhi processes to approximate the Poisson 
jump process. Discretised sample paths can be generated through 
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200 

log Sn+1 = log Sn + i^i-\K- a^/2)5 + aV5Z + ^ NiJi, 



1=1 



where Z denotes a standard normal random variable, Ji ~ N{^j,aj) and 
Ni ~ Bernoulli ((A5/200) exp{-A(5/200}). 

3.3 Double-Exponential Jump Model 

? proposed a jump-diffusion similar Merton's, where the jump size is double-exponentially dis- 
tributed. The double-exponential jump diffusion (DEJD) model is designed to capture the lep- 
tokurtic feature of the empricial return distributions as well as the volatility smile in option mar- 
kets which cannot be successfully modeled by BS-MJ model. The canonical decomposition of the 
driving process of Kou's model is 

(Nit) \ 

dSt = fiStdt + aStdWt + StdiY, -1)1' (12) 

where Wt is a standard Brownian motion, N{t) is a Poisson process with rate A and {Vi} is 
a sequence of independent identically distributed non-negative random variables such that Y = 
log(y) has an asymmetric double exponential distribution with the density 

/y(y) =pr?ie-''^n|,>o} + (1 - p)772e''^^l{,<o}, 
with rji,rj2 > 0, where p>0 represent the probabilities of upward and downward jumps, i.e. 



log(y) = Y=d{ 



i^"*" , with probability p 
— ^~ , with probability 1-p, 



where and ^ are exponential random variables with means rj^ ^ and ^ respectively. In model 
([12]), all sources of randomness, N{t), W{t) and Y's are assumed to be independent. 



The analytical solution of a call option whose price is deteremined by an underlying asset that is 
driven by DEJD model also incorporates a psychological interpretation of investors. As we can see 
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in (jl2p . this model has six parameters, namely fi, the drift parameter, a, the diffusion volatility, A, 
the Poisson rate, p, the probability of having an upward jump, rji, the rate of an upward exponential 
jump and ri2, the rate of a downward exponential jump. By incorporating option prices observed 
with different strikes and maturities, we can improve the estimation, compared with incorporating 
merely the characteristic function of the model. In addition, the option prices used can also enable 
the estimation of the parameters involved in the utility function. 



Given ()12p . one can write down the dynamic of dlogSt by using Ito's Lemma: 



dlogSt 



(^fi - (it + adWt + ^ 



from which we can derive the characteristic function of log St (see ?) under the risk-neutral prob- 
ability measure without taking the jump risk into account: 



21 



log 5o + r - — 



u -|- iu\ 



P 



1 — p 



A+ — iu X- + iu 



since the Levy density of the jump is 



=pXr]ie ''^'^■l{x.>o} + (1 -J')A??2e''2''l{^<o}. 
The corresponding European call price can be obtained via ? method, which is specified as follows: 

27r J-oo a'^ + a — + i[2a + l)u 

Using the independence between the expo- nential and normal distributions used in the model and 
formulae for the sum of double exponential random variables, ? obtains the probability density 
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function of the return, which can be approximated by the following density function: 



1 — X6 f X — fj,6 
-(f) 



X — — a^r]i5 



ay/6 

2^ 



which can be used to define the Radon-Nikodym derivative to adjust for the difference between a 
risk-free probability measure (in Merton's sense) and the physical measure since 



fnix; r) 



The estimating equaiton derived from the option price is 



= 



(e-'-'^ max{e^ - m, 0} - C(5„, T, r)/Sn) (R; 9) 



3.4 With Jump Risk Premium 

? considered a typical rational expectations economy (?) in which a representative investor has 
the utility function of the special form, as in ?: 



U{c,t) = < 



e-'^*-(f'/a , if < a < 1 
e-«* log(c) , if a = 0, 



(13) 



with [/c(c, t) = —^r-^- The goal of the representative investor is to obtain maxc E[J^ U{c{t), t)dt]. 
In his model, Kou also assumed Et, an endowment process, which is, under the physical measure 
P, specified as follows: 

(14) 



^ = fiidt + aidw}^^ + d 



N(t) 

1=1 



15 



given the endowment process ()14p , the asset price wih have the dynamic of the form 



N{t) 



1=1 



(15) 



where dW^"^^ is a Brownian motion independent of dwj^^^ and Vi = vf. Furthermore, a and k in 
(fTSll are related as follows: 



/i = K + (1 - q) <j ^1 - 2^i(^ - a) + o-icrp 



where = ElV" - 1]. 



It can been shown (see, for example, ?) that, under mild conditions, the rational expectations 
equilibrium price, or the "shadow" price, of the security p{t), must satisfy the Euler equation 



(16) 



where Uc is the partial derivative of U with respect to c. To simplify the model, we assume Et = St, 
i.e. pi = p, ai = a and p = f3 = 1. It follows that, as shown in (10) of ?, the Radon-Nikodym 
derivative between the risk-free measure Q and the physical measure P is given by 



— (i?(j+i)5;e) 



^0-+i)^[7e(5(,+i),, (i + 1)5)) 



,(r-K)5 



g<5[(l-a)/.-|a2(l-a)(2-Q))] ^ijj .^ij,(a-l) _ 



c n a— 1 



Here A is a Poisson process with rate A. The jump sizes {11,1^2, • • •} are independent identically 
distributed random variables such that 1^ = log(Vi). The moment generating function of X{t) := 
log{St/So) can be obtained as 



E 



exp{Gie)t}. 



where G{x) = px + ^xV^ + \{E [e^^] - I). In the case of Merton's normal jump-diffusion model. 



_ 1 f x^a'^ 

G{x) = px+ -x^a^ + Xlpjx + - 1 
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and in the case of double exponential jump-diffusion model 

ni \ ~ , 1 2 2 , x/^ , ,\ 

(jix) = iix-\ — X a + A 1 1 

2 X r]2 + x J 

Under the risk- neutral probability Q, we have 



where ( := E [e^] — 1. In the Merton's model 

while in the double exponential jump-diffusion model 

^ ^ PVi _^ (1 - P)V2 _ ^ 
iji-l r?2 + 1 

? adpated the method in ?, which is based on a change of the order of integration, to price 
European call and pution options via Laplace transforms. The Laplace transform with respect to 
k of C{S,e^ ,r,T) is given by 

/oo 
e-^''C{S,e\r,T)dk 
-oo 

^^^exp{G(^ + l)r}, ^>0. 



This leads to the following estimating equation: 

E [ (e-'-'^ max {e^ - m, 0} -c(m,r, ,5,0)) e4(i-")^-|<^'M(2-"))] e«("-i)" 
Note that the option price c under the double exponential jump diffusion dynamics can also be 
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obtained directly using the method proposed by ?, 

c{m;r,6,0) = T + ^a"^ - XC,a,X,p,fji,fj2;log{l/m),5^ 

+ me'^'^T ^r- - X(,a,X,p,fji,fi2;log{l/m),5^ , 

where the definitions of X,p,fji,fj2 and T(-) can be found in ?. 

The characteristic function of a return drived from the price driven by the process ()15p can be 
obtained similarly as in Sections 3.2 and 3.3. The corresponding estimating equations are thus 

E[cos{tR)] - Re((/>(t; 6)) = E[sm{tR)] - lm{(j){t; 0)) = 0. 

4 Asymptotic Results 

Regularity conditions: 

1. E[g{t, Xi; 6Q)g{t, Xi] Oq)'] is positive definite for t G [—a, a], a > 0; 

2. -^g(t, x; 6) is continuous in a neighbourhood of 6q, for t G [—a, a], x G M; 

3. supege \\-§Qg{t,x;e)\\ < H{t,x), where 




H{t,x)dF{x)dGi{t) < oo; 



4. The rank of E[-^g{t, Xi; 6q)] is min{2, d} for all t G [—a, a], where d is the dimension of 0; 

5. ggggi dit, ^'i ^) IS coutinuous in for G 0, i G [—a, o] and x G M; 

6. sup^gQ \\^^)gTg{t,x;6)\\ < H{t,x), where H is given in[3l 

Proposition 4.1. Under conditions 1-4, with probability one, denote 9 = argmin0Ti(0) which 
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satisfies \\G — Oo\\ < n 



Qin{t,e,Xi{t;e)) = 
Q2nit,e,Xi{t;e))dGi{t) = 0, 



where 



g{t;xj,e) 



n ^1 + X'git;xj;i 



Proposition 4.2. Under conditions 1-6, for the estimator 9 given in Proposition 2.1, we have as 
n — >■ oo, 



V^{e - Go) 



S2l{t)s^l{t)Si2{t)dGi{t) 



-1 



S2i(i)«n (*)V^Qin(i;6>o,0)dGi(t) \ + Op{l) 



where 



Suit) 

Sl2{t) 



-E[g{t,Ri;eo)g{t,Ri;9oy], 
s2iit) = E[^9{t,Ri;eo)], 



S2l{t)s^l{t)si2{t)dGi{t) 

ra pa 

/ / S2l{ti)s^^iti)T{ti,t2)s^i\t2)si2it2)dGl{ti)dG2{t2) 
J —a J —a 

X< / S2l{t)s^i\t)si2{t)dGi{t) 



where 



T{ti,t2) = E[g{ti,Rt; eo)g{ti,Rt] 0o)']- 



Corollary 4.1. When ki > k2, the asymptotic variance S = of \/n(6 — 6) cannot decrease if 
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an estimating equation is dropped. 



5 Numerical Results 
5.1 Simulations 

For each model, 500 sample paths with size n = 125, 250, 500 and 1000 starting at initial value 
log 5*0 = 100 with frequency 6 = 1/52 were simulated. Similar to Chan et al. (2009) approach, 
we also choose the uniform weight function G{t) and In can be approximated by the Riemann sum 
of ln{t) evaluated at t G [—5.0,5.0] with the number of grids set to be 100. |1| Simulation results 
for BS, BS-MJ and DEJD are tabulated in Tables [T|2] and HI respectively. As we can see from the 
simulation results, by incorporating more option prices, the estimated standard deviation of the 
estimates are reduced, which is due to the result of Corollary 14. 1[ 



5.2 Case Study 

We examine empirically whether the proposed methodology can be applied to the real data set 
and what insights call prices can reveal when we incorporate them into the model. Historical S&P 
500 index values and corresponding call option prices between 2 Janurary 1987 and 31 December 
2008 were downloaded from Wharton Research Data Services (WRDS). The index is sampled at 
4-day frequency and in total we have 1,260 data points. The duration of four days is chosen so as 
to match the time to maturity of the call prices. The calls were traded on Chicago Board Options 
Exchange (CBOE). 

We included, in our simulations, from one to four call prices that were most frequently traded 
on the last day of our analyses so as to reflect the market information on that particular trading 
day. The mean annual rate of return is 0.0531 with the associated volatility equals 0.1328. In 
addition to the market crash of 1987, the tech- and credit-bubble between the late 90's and mid 
2000's as well as September 11 attack in 2001, the sample period also covered the recent Lehmann 

^In Chan et al. (2009), they chose the interval to be [—0.5,0.5]. In the simulation studies, we found that using 
[—0.5, 0.5] produced poor estimations. A wider interval chosen allows the data to provide more information about 
the parameter values. 
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Brother's collapse as a result of credit crunch in 2008. In particular, our data analysis was done 
with the last day selected as September 29, 2008 - the day on which the largest single day plunge 
was recorded shortly after Lehman brothers' and Washington Mutual's bankrupcy. Furthermore, 
on that day, the Volatility S&P (VIX), a measure of market volatility, has the record highest jump 
in history. Estimated values of the parameters and the associated esimated asymptotic variances 
are tabulated in Tables 6-8. 

It can be seen from the tables that by incorporaing constraints due to observed option prices, 
one can lower the variance of the estimates. It should be also noted that in order for ?'s approach to 
achieve the same magnitude of variance as what we can see by including additionally one call price, 
the sample size should have to be roughly doubled. In other words, using call prices as constraints 
reduces the required sample size at the expense that the equity price dynamics arc specified by a 
particular model. Since the option prices are considered as a summary of the current market view 
on the underlying equity price dynamics, our methodology can successfully capture more updated 
estimate of the current market condition. This can be seen in the data analysis results in which the 
volatility and/or jump size estimates are both larger than the estimates that ? provided, which can 
be interpreted as the consequence because of the late-2000's financial crisis. Finally, we comment 
that, due to the small number of option prices included, it is challenging to produce an accurate 
estimate a, the risk-preference parameter of investors of which informaiton can be only derived 
from the option prices. 

6 Conclusion 

Levy processes arc an excellent tool for modeling price processes in mathematical finance. Its pop- 
ularity arises from its flexibility and simple structure in comparison with general semimartingales. 
Estimation for Ley processes are challenging statistical inference problems because of the lack of 
analytical expression for the transitional density function. Inspired by ? that uses integrated em- 
pirical likelihood approach for parameter estimation, we propose in this paper incorporating call 
prices as constraints in addition to using the characteristic function associated with the process. 
This method provides a more efficient estimate that can reflect the recent market condition more 
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accurately which is demonstrated via simulations and real data analyses. The idea of using deriva- 
tive prices as one of the estimating equations is not restricted to call prices only; in fact, any price 
that can be expressed in terms of expectation of an independent random variable that follows the 
same distribution as specified by the underlying process are eligible for being included as one of 
constraints. The approach, therefore, has robust theoretical and versatility for a wide range of 
processes including processes with jump components. 

7 Appendix 



Proof of Proposition 4-1- It follows closely the proof of lemma 1 of ?. 



□ 



Proof of Proposition 4-^- Similar to ?, we can show that 



^Qin{t;e;0) ^ Si2{t), 

-^Qin{t;e;0) 4 sn{t), 

^Q2n(t;^;0) = 0, 

^Q2n{t;e;0) 4 S2l{t) 



uniformly in t E [—a, a]. Denote \\6 — Oo\\ + supjg[_(j „] ||Ai(t;0)||. Then, we can expand Qin and 
Q2n using Taylor series expansions and yield 



WQin(i;6>o,0) - s^l{t)si2{t){9 - Oo) + Op{6n), and 
Q2n{t;e,\iit;e))dGiit) 

Q2nit; 00, 0) + ^9^ll(^^^0 _ Oo) + ^^4^A,(t; 9) \ dG,{t) + o,(5. 



de 



dX' 



J S2i{t)s^^\t)si2{t)dGiit)'^ xjy" S2iit)s^^\t)Qinit;Qinit;eo,0))dGi{t) 

+Opi5n), 



which completes the proof. 



□ 
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Proof of Corollary 5.1. To prove the Corollary, it suffices to show 

/ / S2i{ti)s^^{ti)T{ti,t2)s^i\t2)si2{t2)dGi{ti)dGi{t2)- [ S2i{h)s^l{h)si2{ti)dGi{ti) < 0, 

J —a J —a J —a 

where A < B denotes A—B isa, negative-semidefinite. Observe that, by Cauchy-Schwarz inequahty, 



pa pa 

/ / S2i(ii)sn(ii)r(ii,i2)sn (*2)si2(i2)dG'i(ii)dGi(i2) 

J —a J —a 



< 



hi(ti)sn (ti)sii(^i)sii(ii)s2i(ti)]^^^ [sii(t2)sn (t2)si2(^2)sii(^2)s2i(i2)]'^'dG'i(ti)o!Gi(t2) 



1/2 



= r (s2i(ti)sri'(ii)s2i(ti))'^'dGi(ti) 

y —a 

< [ S2l{h)s^l{h)s2l{tl)dGi{ti). 
J— a 



It follows that E < 
Write 



J^^S2i{t)s^l{t)s2i{t)dGi{t) . The proof can be completed following ?. 



si2{t;e) = 
sii{t;e) = 



dgiV f dgk- 

ae J ' ■ ■ ■ ' V de 

Sll,ait,0) Su^h{t,0) 

sii,c{t,0) sn,d{t,9) 



dgk 
06 



where sii,a(i, 0) is a (fc — 1) x (fc — 1) matrix. Then, for all t G [—a, a], 



S2i(i)sn(i)s2i(t) = 



Si2(i; 9), 



> 



Si2{t;9), 



dgk^ 






de J 


)] 






. / 


dgk^ 






de , 








V 



sii,a{t,e) sii^b{t,e) 
sii,c{t,e) su^d{t,e) 



si2it;0) 



dgk 

de 



Sll,a(t,6') 



= s^2{t;e)sn,a{t,0)s^2{t;ey 



dgk 

de 



which completes the proof. 



□ 
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strike 


1 strike 


2 strikes 


4 strikes 


K 


NA 

(Chan et al. 2009) 


0.995 


0.995, 1.015, 


0.98, 0.995 
1.015, 1.025 


n = 125 


fi = 0.044(0.192) 

d = 0.298(0.019) 


l2 = 0.056(0.190) 
a = 0.2992(0.012) 


ju = 0.045(0.189) 
a = 0.2998(0.009) 


/I = 0.041(0.117) 

d = 0.299(0.008) 


n = 250 


J2 = 0.047(0.133) 
a = 0.299(0.014) 


/i = 0.054(0.132) 
a = 0.2996(0.009) 


/I = 0.050(0.130) 
d = 0.2999(0.0088) 


= 0.042(0.130) 
a = 0.2869(0.073) 


n = 500 


l2 = 0.051(0.097) 

a = 0.300(0.010) 


/I = 0.052(0.091) 
a = 0.2998(0.0070) 


/2 = 0.051(0.0984) 
a = 0.2999(0.0071) 


l2 = 0.053(0.0873) 

a = 0.2996(0.0064) 


n = 1000 


J2 = 0.047(0.068) 
a = 0.3000(0.007) 


/2 = 0.054(0.069) 
a = 0.2949(0.0069) 


fi = 0.054(0.0069) 
a = 0.295(0.0066) 


/I = 0.052(0.0071) 
d = 0.295(0.0062) 



Table 1: Black-Scholes Model (BS) with true values (/x,cr) = (0.050,0.30). 





strike 


1 strike 


2 strikes 


4 strikes 


K 


NA 

(Chan et al. 2009) 


0.995 


0.995, 1.015 


0.98, 0.995 
1.015, 1.025 


n= 125 


^ = 0.068(0.5829) 
a = 0.1597(0.3086) 
A = 2.7485(10.0591) 
■jlj =-0.3121(0.6929) 
^} = 0.4440(0.4946) 


/2 = 0.0329(0.1903) 
a = 0.2305(0.1608) 
A = 2.8674(6.0220) 

/I} =-0.3445(0.4306) 
aj = 0.4665(0.2547) 


/2 = 0.0103(0.1281) 
a = 0.2392(0.1602) 
A = 3.1887(5.4130) 
■fl} =-0.2749(0.2912) 
(j} = 0.4474(0.2468) 


/2 = 0.0244(0.1644) 
a = 0.2709(0.1378) 
A = 3.0066(2.8460) 
flj =-0.2670(0.2848) 
= 0.5009(0.2273) 


n = 250 


/2 = 0.0639(0.3790) 
a = 0.1613(0.2518) 
A = 1.8401(3.5736) 

fTj =-0.3172(0.4775) 
aj = 0.4229(0.4559) 


/i = 0.0443(0.1435) 
a = 0.2304(0.1424) 
A = 2.3803(2.9758) 
JIJ =-0.2993(0.3115) 
aj =0.5126(0.2216) 


/2 = 0.0258(0.085) 
a = 0.2561(0.1303) 
A = 2.5194(2.0394) 

JI} =-0.2811(0.2445) 
a} = 0.5047(0.2254) 


/2 = 0.0235(0.0824) 
a = 0.2806(0.1250) 
A = 2.5063(2.0116) 

jTj =-0.2723(0.2397) 
aJ = 0.5269(0.2057) 


n = 500 


/i = 0.0740(0.2405) 
a = 0.2033(0.2094) 
A = 1.7045(0.9860) 
JIJ =-0.2411(0.3034) 
a} = 0.5153(0.2159) 


/i = 0.0650(0.1228) 
a = 0.2259(0.1379) 

A = 2.0594(1.4227) 
pTj =-0.2864(0.2636) 
a} = 0.5203(0.2143) 


/I = 0.0524(0.1057) 
a = 0.3290(0.1097) 

A = 1.7884(1.2292) 
fTj =-0.3097(0.2377) 
a} = 0.5344(0.2119) 


/I = 0.0597(0.1076) 
d = 0.3498(0.1226) 
A = 1.9321(1.0289) 
JJ} =-0.3049(0.2280) 
a} = 0.5748(0.1987) 


n = 1000 


fi = 0.0315(0.1900) 
a = 0.2349(0.1658) 
A = 1.8302(0.9038) 
fj =-0.2807(0.2423) 
aj = 0.5800(0.2155) 


l2 = 0.0317(0.0884) 
a = 0.2553(0.1282) 

A = 1.9637(0.8853) 
jTj =-0.3068(0.2161) 
aj = 0.6030(0.1720) 


ju = 0.0326(0.0825) 
a = 0.3432(0.1017) 
A = 1.6987(0.7044) 
=-0.2999(0.1857) 
aj = 0.6095(0.1618) 


l2 = 0.0396(0.0813) 
a = 0.3407(0.0975) 
A = 1.9238(0.7021) 
JjTj = -0.2776(0.1627) 
aJ = 0.6051(0.1470) 



Table 2: Black-Scholcs model with Merton Jumps (BSMJ) with true values {fj,,a,X,fj,j,aj) = 
(0.0875,0.30,2.0,-0.2,0.60). Notice that in this simulation study, we do not follow Chan et al. 
(2009) true parameter values because it is very unlikely to have jump size with mean and variance 
2.0 and 6.0 respectively while the underlying drift and variance are mere 0.05(5 and 0.30^. 
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1 strike 


2 strikes 


4 strikes 


K 


0.995 


0.995 
1.015 


0.985,0.995 
1.015,1.025 


n = 125 


/I = 0.0974(0.1840) 
a = 0.2607(0.1305) 
A = 1.9939(0.0349) 
JJT} = -0.0874(0.1070) 
aj = 0.0370(0.1530) 
a = 0.6956(1.4105) 


/I = 0.0794(0.1682) 
a = 0.2741(0.1193) 
A = 1.9956(0.0306) 
■jTj = -0.1004(0.0994) 
ct} = 0.0564(0.1612) 
a = 0.6460(1.2068) 


/I = 0.0586(0.1603) 
a = 0.2826(0.1069) 
A = 1.9980(0.0311) 
JJ} = -0.1088(0.1051) 
5v = 0.0592(0.1642) 
a = 0.6429(1.1079) 


n = 250 


= 0.0966(0.1422) 
a = 0.2792(0.1140) 
A = 1.9934(0.0250) 
jTj = -0.0811(0.0847) 
a} = 0.0461(0.1668) 
a = 0.6744(1.1766) 


= 0.0866(0.1382) 
a = 0.2789(0.1053) 
A = 1.9920(0.0511) 
jTj = -0.0927(0.0818) 
a} = 0.0496(0.1659) 
a = 0.6382(1.0552) 


/2 = 0.0732(0.1315) 
a = 0.2921(0.0756) 
A = 1.9955(0.0247) 
/i} = -0.1009(0.0853) 
a} = 0.0523(0.1641) 
a = 0.6449(0.9589) 


n = 500 


/2 = 0.1027(0.1121) 
d = 0.2949(0.0699) 
A = 1.9938(0.0203) 
fij = -0.0841(0.0693) 
aj = 0.0607(0.1655) 
S = 0.5703(0.7906) 


= 0.0995(0.1087) 
d = 0.2982(0.0552) 
A = 1.9932(0.0154) 
fj = -0.0897(0.0644) 
a} = 0.0605(0.1635) 
S = 0.5499(0.7051) 


/2 = 0.0876(0.1047) 
d = 0.2982(0.0391) 
A = 1.9939(0.0136) 
iTj = -0.0990(0.0623) 
3v = 0.0819(0.1608) 
a = 0.5846(0.6543) 


n = 1000 


/2 = 0.1044(0.0860) 
a = 0.3032(0.0461) 
A = 1.9946(0.0096) 
fTj = -0.0822(0.0594) 
a} = 0.1007(0.1488) 
a = 0.5846(0.6543) 


/I = 0.1040(0.0869) 
a = 0.3026(0.0369) 
A = 1.9947(0.0080) 
■jTj = -0.0832(0.0551) 
a} = 0.1042(0.1493) 
a = 0.5678(0.6310) 


/2 = 0.0918(0.0855) 
a = 0.3040(0.0.0279) 

A = 1.9960(0.0068) 
fi} = -0.0946(0.0494) 
aj = 0.1212(0.1386) 

a = 0.5722(0.5827) 



Table 3: Merton Jump-diffusion model with true values {ii,a, X, iij,aj,a) 
(0.095, 0.30, 2.0, -0.08, 0.20, 0.60). 
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1 strike 


2 strikes 


4 strikes 


K 


0.995 


0.995 
1.015 


0.985,0.995 
1.015,1.025 


n = 125 


Ji = 0.1592(0.1997) 
a = 0.2274(0.1999) 
A = 1.7825(0.9630) 
p = 0.1107(0.7369) 
rfi = 7.6935(0.3580) 
m = 8.8874(0.7513) 
a = 0.3415(1.7920) 


ju = 0.1473(0.1713) 

a = 0.3109(0.0658) 
A = 1.8095(0.6815) 
p = 0.2571(0.9921) 
fji = 7.6935(0.3580) 
rfs = 8.9638(0.5798) 
a = 0.1705(1.3924) 


fi = 0.1129(0.1777) 
a = 0.3141(0.0528) 

A = 1.9343(0.6037) 
p = 0.3422(0.5358) 
rfi = 7.5872(0.4652) 
f]2 = 8.9658(0.4665) 
a = 0.2717(1.2334) 


n = 250 


/2 = 0.1152(0.1876) 
a = 0.2969(0.0641) 
A = 1.8130(0.6496) 

p = 0.3424(0.4645) 
rfi = 7.5572(0.6307) 
r]2 = 8.9710(0.7906) 
a = 0.1835(1.2546) 


/I = 0.1115(0.1527) 
a = 0.3018(0.0443) 
A = 1.8397(0.6693) 

p = 0.3934(0.4164) 
rfi = 7.5784(0.3063) 
% = 9.0307(0.4292) 
a = 0.0800(1.1524) 


J2 = 0.1007(0.1469) 
a = 0.3023(0.0606) 
A = 1.9636(0.5531) 
p = 0.3897(0.3653) 
?fi = 7.5292(0.4211) 
f]2 = 8.9892(0.4004) 
a = 0.2256(1.0019) 


n = 500 


/2 = 0.1044(0.1478) 
a = 0.3065(0.0355) 
A = 1.8205(0.6064) 
p = 0.4018(0.3516) 
ffi = 7.5446(0.3409) 
rf2 = 9.0713(0.3104) 
a = 0.0764(0.8658) 


/} = 0.1206(0.1239) 
d = 0.3036(0.0329) 
A = 1.8907(0.7034) 
p = 0.4041(0.3138) 
rfi = 7.4794(0.4696) 
rf2 = 9.0671(0.3318) 
a = 0.0220(0.8857) 


/i = 0.0929(0.1245) 
a = 0.2976(0.0487) 
A = 2.0139(0.5391) 
p = 0.4030(0.3050) 
?fi = 7.4912(0.3298) 
% = 8.9822(0.2867) 
a = 0.1735(0.7628) 


n = 1000 


= 0.1109(0.1201) 
a = 0.2889(0.0366) 

A = 2.0648(0.6382) 
p = 0.4004(0.2498) 
rfi = 7.4740(0.2350) 
r]~2 = 8.9613(0.3358) 
a = 0.1242(0.6115) 


/2 = 0.1169(0.1031) 
a = 0.2891(0.0400) 
A = 2.1036(0.7695) 
p = 0.4026(0.2132) 
rfi = 7.4330(0.3336) 
ff2 = 8.9589(0.3884) 
a = 0.0254(0.5691) 


/2 = 0.0725(0.0955) 
a = 0.2982(0.0343) 

A = 2.0227(0.4550) 
p = 0.4398(0.2069) 
jfi = 7.4613(0.3112) 
r}~2 = 9.0012(0.2128) 
a = 0.1828(0.4250) 



Table 4: Kou Double-exponential Jump-difFusion model with true values (/i, cr, A,p, r/i, 772, a) = 
(0.095, 0.30, 2.0, 0.05, 7.5, 9.0, 0.20). 



strike 
Chan et al. (2009) 


1 strike 


2 strikes 


4 strikes 


/2 = 0.0620(0.0296) 
a = 0.1327(0.0036) 


/2= -3.5483 X 10-^'J(0.0179) 
a = 0.2058(0.0025) 


/2 = 0.0133(0.0030) 
a = 0.2257(0.0019) 


fl = 0.0133(7.7640 x lO^*^) 
a = 0.2340(6.2043 x 10"^) 



Table 5: Empirical estimation for the S&;P500 index between January 2, 1987 and September 29, 
2008: Black-Scholes Model. 
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strike 
Chan et al. (2009) 


1 strike 


2 strikes 


4 strikes 


/2 = 0.0800(0.0296) 
a = 0.12715(0.0137) 
A = 1.24319(0.0049) 
-fTj = 0.0213(2.3519) 
aj = 0.0267(2.5951) 
NA 


/I = -0.0763(0.0031) 
a = 0.2101(0.0177) 
A = 1.8688(0.0002) 

fj = 0.0177(0.9440) 
aj = 0.0821(0.943) 
a = 0.2860(2.2553) 


/i = -0.0731(0.0009) 
a = 0.2374(0.0080) 
A = 1.8760(0.0011) 
fj = -0.0050(0.0622) 

3v = 0.0223(0.2606) 
a = 0.2889(1.5988) 


fl = -0.0755(0.0005) 
a = 0.2379(0.0048) 
A = 1.876(0.0001) 
JIJ = -0.0061(0.0283) 

aj = 0.0226(0.1270) 
a = 0.2934(1.3815) 



Table 6: Empirical estimation for the S&P500 index between January 2, 1987 and September 29, 
2008: Black-Scholes with Merton Jumps Model. 



strike 
Chan et al. (2009) 


1 strike 


2 strikes 


4 strikes 


= 0.0606(0.0279) 
a = 0.0924(0.0048) 
A = 2.160(0.1067) 
p = 0.2715(0.0001) 
?fi = 16.1389(0.2696) 
r]2 = 25.2204(1.6940) 
NA 


ll = -0.0438(0.0173) 
a = 0.09447(0.0043) 
A = 2.1493(0.0828) 
p = 0.2408(0.0001) 
jfi = 16.1391(0.2761) 
r]2 = 25.2201(1.5885) 
a = 0.09093(7.9389) 


/2 = 0.0035(0.0034) 
a = 0.1956(0.0014) 
A = 2.1684(0.0371) 
p = 0.4576(0.0001) 
rfi = 16.1434(0.0775) 
rf2 = 20.1736(0.5765) 
a = 0.2084(0.4855) 


J], = -0.0775(0.0014) 
a = 0.2246(0.0014) 
A = 1.6002(0.0323) 
p = 0.4633(0.0001) 
rfi = 21.5180(0.0589) 
rf2 = 30.2644(1.3122) 
a = 0.2329(0.1077) 



Table 7: Empirical estimation for the S&:P500 index between January 2, 1987 and September 29, 
2008: Kou Double Exponential Jump-Diffusion Model. 
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